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Q ■ Abstract 

■ We compute the leading order contribution to the stress-energy tensor corresponding to the 

, modes of a quantum scalar field propagating in a Priedmann-Robertson- Walker universe with 

I arbitrary coupling to the scalar curvature, whose exact mode functions can be expanded as an 



O 



infinite adiabatic series. While for a massive field this is a good approximation for all modes 
when the mass of the field m is larger than the Hubble parameter H, for a massless field only 



the subhorizon modes with comoving wave-numbers larger than some fixed obeying fc^/a > H 
can be analyzed in this way. As infinities coming from adiabatic zero, second and fourth order 

> 

■ expressions are removed by adiabatic regularization, the leading order finite contribution to the 

, stress-energy tensor is given by the adiabatic order six terms, which we determine explicitly. For 

in 

T^lj- I massive and massless modes these have the magnitudes H^/m^ and H^a^/kl, respectively, and 



higher order corrections are suppressed by additional powers of {H/rnf' and (Ha/k^:)'^. When the 
scale factor in the conformal time r/ is a simple power a(ry) = (1/r/)'^, the stress-energy tensor obeys 



^ I P = ijjp with uj = [n — 2)/n for massive and w = (n — 6)/(3n) for massless modes. In that case, 

the adiabaticity is eventually lost when < n < 1 for massive and when < n < 3/2 for massless 
fields since in time H/m and Ha/k^ become order one. We discuss the implications of these results 
for de Sitter and other cosmologically relevant spaces. 
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I. INTRODUCTION 



Determining vacuum energy in quantum field theory for a given physical situation is an 
important problem which may lead one to deduce significant theoretical and observational 
results. In the simplest case where one considers a free quantum field confined in between 
two parallel plates in fiat space, the existence of Casimir energy is experimentally verified 
and it agrees with the field theory calculations. This motivates one to search for possible 
cosmological impacts of vacuum energy at early or late times. Especially with the discovery 
of the recent accelerated expansion of the universe, it is important to see whether vacuum 
energy can play a role in acceleration. Indeed, the cosmological constant is usually thought 
to be related to the vacuum energy (see e.g. [1]). Of course, the problem of fixing vacuum 
energy in a cosmological setting is more complicated than determining the Casimir energy 
associated with parallel plates. Conceptually, the most important difference is the non- 
uniqueness of the vacuum in an expanding universe due to absence of Poincare symmetry. 
Regularization in a curved space-time is also more subtle and difficult compared to the the 
flat space examples since it depends on the geometry in a non-trivial way. 

Adiabatic regularization [2-4] is a convenient way of obtaining finite results out of diver- 
gent stress-energy tensor expressions in Friedman-Robertson- Walker (FRW) space-times. In 
that scheme, one subtracts mode by mode contributions to the stress-energy tensor coming 
from a suitably defined adiabatic basis. This is in principle similar to subtracting infinite 
flat space contribution to get a flnite Casimir energy for parallel plates. One nice feature 
of adiabatic regularization is that the final stress-energy tensor is guaranteed to be con- 
served. Moreover, it is known to be equivalent to point-splitting regularization in FRW 
space-times [5, 6]. Adiabatic regularization may suffer from infrared divergences for mass- 
less fields [7], so extra care is needed in such cases. However, it is possible to obtain, for 
example, the standard trace anomaly in a relatively simple way [8] , which gives further con- 
fidence to the method. To remove quartic, quadratic and logarithmic ultraviolet divergences 
that generically appear in the stress-energy tensor, it is enough to determine the adiabatic 
mode functions up to fourth order time derivatives (i.e. up to adiabatic order four) and the 
corresponding subtraction terms are determined in [9]. 

To specify the vacuum state in a FRW space-time, one should fix the time dependence 
of the mode functions in a certain way. The most convenient choice is the so called Bunch- 
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Davies vacuum, where one identifies the "negative frequency" solution in the remote past 
with the mode function corresponding to the annihilation operator (there is however an in- 
herent ambiguity in determining the vacuum if one considers a realistic cosmological scenario 
[10].) As we will discuss in the next section, sometimes the mode functions can be solved 
as a well defined infinite series in the adiabatic expansion scheme, i.e. instead of stopping 
at order four to get an approximate adiabatic solution, which would give the adiabatic sub- 
traction terms for regularization, one can in principle continue to obtain an exact solution 
in series form. The vacuum associated with such mode functions is called the adiabatic 
vacuum} Naturally, the corresponding stress-energy tensor can be calculated as a series 
and now adiabatic regularization requires throwing out the adiabatic zero, second and the 
fourth order contributions, leaving the sixth order terms as the leading order contribution 
to the finite stress-energy tensor. Due to presence of adiabaticity, the stress-energy tensor 
can be thought to be related to vacuum polarization effects rather than the particle creation 
ones. 

In this paper, we explicitly calculate adiabatic order six terms for the stress-energy tensor 
of a scalar field propagating in a FRW space-time with arbitrary coupling to the curvature 
scalar. As we will discuss, this is a good approximation for a massive field if the mass 
is larger than the Hubble parameter. For a massless field only the contributions of modes 
which have sufficiently large (comoving) wavenumbers can be determined in this way. As one 
would expect, the sixth order expressions are complicated and thus not very illuminating. 
However, when the scale factor of the universe is a simple power in conformal time, the 
stress-energy tensor simphfies considerably. We elaborate on possible implications of our 
results for cosmology. Specifically, we fix the magnitude (and the sign) of the vacuum 
energy density corresponding to adiabatic vacuum in cosmologically relevant spaces, and 
determine when the assumption of adibaticity is a suitable approximation. 



II. ADIABATIC VACUUM 



We consider a scalar field (f) which has the following action 



/^[{V4>f + (m^ + iR)4>'']- (1) 

^ One usually defines adiabatic vacuum up to a certain order by identifying the initial values of the exact 
mode functions with the approximate adiabatic mode functions determined to that order. The adiabatic 



vacuum we define here has infinite order. 



The coupling of the scalar field to the curvature scalar is governed by the dimensionless 
parameter ^. Varying the action with respect to the metric one can determine the stress- 
energy-momentum as 

= V^0V,0 - eV^V,(02) - ^g^, [(V0)2 + (m^ + CR)<f>' - 2ev2(0')] + e^'i^M- (2) 

As our background we take the FRW space-time with the metric 

ds'' ^a{r)f{-dr)'' + dx^ + dy^ + dz''). (3) 



For later use we define the Hubble parameter in conformal time as 



a 



(4) 



where the prime denotes derivative with respect to 77. The quantization of the scalar field 
in a FRW background is straightforward. Defining a new field fi by 



IJ, = a( 



(5) 



and applying the standard canonical quantization procedure, one can see that the field 

operator n can be decomposed in terms of the Ume-mdependent ladder operators as 

d^k 



(27r)3/2 



/i,(77)e^'=-a,- + //,(77)*e-'=-al 



(6) 



where k is the comoving momentum variable, [a^, at ] — 5{k — k') and the mode functions 
satisfy the Wronskian condition ^lk^Jl'k ~ l^ll^k — ^ together with 



1^1 + 



r + m^a^ + (6^- 1)- 



l^k = 0. 



The ground state |0 > of the system can be defined by imposing 



agio >= 0. 



(7) 



(8) 



Using the expression for the stress-energy- momentum tensor (2), one can calculate the vac- 
uum expectation values as^ 



1 r°° r 9 / \ 

< 0|p|0 > = [\iJi'^ - hiJikV + [k'' + m'a' - 6^/1') + 6^/1(1//^!')' 



k'^dk, (9) 



< OIPIO > 



_ I r 



li'j^-h/ikl - { — + m a ] \lJ'k\ +6^/i(|/Xfe 



2\/ 



2\ll 



k'^dk. 



As it is written in (2), there is no ordering ambiguity in the stress-energy- momentum tensor operator. 
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The system is now fully specified except the mode function /i^ obeying (7), which would 
also define the vacuum state |0 >. 

For the adiabatic expansion, one writes /ik as 

Written in this way, the mode function /i^ automatically satisfies the Wronskian condition. 
Using (7), Qk can be seen to obey 



a ^ 4nl 2 Ofe 



k' + m'a' +(6^-1)- + ^7^-^^. (11) 



It is possible to solve the above equation iteratively as follows: One starts from the zeroth 
order solution that contains no time derivatives, i.e — \^^W+rnFa?. Using in the 
right hand side of (11), one can determine a second order solution fij,^', which contains terms 
up to two time derivatives. Now Qf, can be used in the right hand side and this procedure 
can be continued iteratively to get a series solution for ^2^. It is important to emphasize that 
in this series expansion the number of time derivatives acts like a perturbation parameter. 
Assuming that the final infinite sum converges, one gets a unique function /x^. To obtain 
the second linearly independent solution, one should start the series with the negative root 

Let us try to see when the above prescription can give a well defined and thus a 
solution for n^. For that let us analyze the second order solution which can be found as 

((3(-l)«" m\(,'^ + aa") 5ni^(i^a'^ 



(12) 



2a{k^ + m^a'') 4:{k^ + m^a^)^ 8{k^ + m^a^)\ 

For a massive field, one sees that the second order terms in the square brackets have their 
largest vahics for A; = 0, which have the magnitude H'^/m^, where H is the Hubble parameter 
with respect to the proper time 

H--,. (13) 

Therefore, as long as m ^ H, the second order terms will be much smaller than the zeroth 
order ones even for the k — mode (note that the corrections are more suppressed for larger 
k). By inspecting the higher order adiabatic contributions, one can see that the magnitude 
of the n'th order adiabatic terms is equal to (i?/m)" for the k — mode. Thus, lor H 
the adiabatic expansion is trustable for all modes to determine /x^. 



On the other hand, for a massless field with m — 0, the second order solution becomes 



[2] 



k 



1 + 



(6e-l)a" 
2aP 



(14) 



This time higher order adiabatic corrections are suppressed by powers of Ha/k. Not sur- 
prisingly, adiabatic expansion fails for modes with k/a < H (i.e. for superhorizon modes). 
However, the expansion can still be used to determine the mode functions for k/a ^ H 
(i.e. for subhorizon modes). Since the vacuum is defined mode by mode for each fc, one can 
use adiabatic expansion to fix with large enough comoving wavenumbers with k > k^ for 
some fixed k^ obeying 

K > aH. (15) 

Moreover, the momentum integrals in (9) can be decomposed into two decoupled pieces 
corresponding to the intervals (0, /c*) and (fc^,, oo), where in the second interval the adiabatic 
expansion can safely be used. 

We thus conclude that the adiabatic vacuum is physically viable^ for all modes of a 
massive scalar field \lm^ H and it can only be imposed for the ultraviolet (UV) modes of 
a massless scalar obeying A; > A;* ^ aH, for some fixed k^. 

Since the mode functions are (uniquely) specified by the adiabatic expansion scheme, 
one can calculate the vacuum expectation values (9) for the adiabatic vacuum. For this 
calculation, it is convenient to express (9) in terms of Vt^. Using (10) one finds 



< 0|p|0 > 



< OIPIO > 



1 



1 

4^ 



VLk 1 

— + 

2 20^ 

fifc 1 

— + 

2 20j^; 



h + 



h + 



20t 



+ 



^^{e + m'a'-6Ch')-3^h^ 



k^dk, 



20. 



nf 



k'^dk. (16) 



To obtain the leading order contribution to the stress-energy tensor one should use the sixth 
order solution = and furthermore subtract the infinite zeroth, second and fourth 
order adiabatic terms which can be determined by using 0,^ — a-g^'in in (16). It is 
important to recall that in this whole procedure the number of time derivatives acts like a 
perturbation parameter. 

^ In general one may be concerned with the existence of rij.°°' since one actually makes an asymptotic 
expansion about a non-analytical point and the convergence of the series may fail in a very short time 



[10] . This should not be an issue for very massive fields or for modes with very large wave-number. 
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This is a straightforward but a very cumbersome calculation to carry out, which we 
perform with the help of a computer. The final result for Jl'^l is very comphcated. However, 
after using Q'^l in (16) and performing the elementary and convergent momentum integrals, 
we obtain for the massive field the following relatively simple expression^ 

p{m) = 40320^,2^2^12 [(844 - 35280a'^ - 72(173 - 1596^ + 3780^') aa'^a"+ 
12(683 - 6636^ + 16380^^)a2a'=^a'" + 9(1237 - 14924^ + 59220^^ - 75600^3)0^0'^ a"^ 
+9(17 - 168e + 420e')(-16a=^a'2a(^) + o^o^^^' - 2o^o"o(^) + 2o^o'o(^)) 
+(673 - 8316^ + 34020^2 - 45360^^) (2o=^o"=^ - 6o=^o'o"o(3))] , 

P{m) = 40320^,2^2^12 [12(211 - 8820a'' - 8(4363 - 39186^ + 90720e')aa'^o" 

+4(5819 - 56028^ + 137340^^) o^o'^o^^^ + 15(2389 - 31108^ + 107100^^ - lQ584Q^^)a^a'^a"^ 

+10(-473 + 5628^ - 21924^^ + 27216^3)o=^o"=^ + 4(617 - 6930^ + 23940^^ - 22680^=^) o^o"o(^) 

+ (1601 - 19152^ + 74340^2 - 90720^3)a^a(3)' - 6(17 - 168^ + 420^2)a^a(^) 

+6(-3949 + 44828^ - 160020^^ + 166320^3)o3o'o"o(3) + 4(1907 - 18732^ + 46620^2)o3o'^o(^) 



+78(17 - 168^ + 420^^)o^o'o(^) 



(17) 



where the numbers in the parentheses which appear above the scale factor indicate the order 
of r) derivatives. 

As discussed above, for a massless field only modes with k > K can be treated adiabat- 
ically. In that case, one can calculate the partial contribution of these modes to the total 
stress-energy tensor by changing the limits of the /c- integrals in (9) to the range {k^,oo). 
Note that this partial stress-energy tensor totally decouples from the rest of the modes and 
it is self consistently conserved. For the massless case we then find 
(1 - 60' 



256A;27r2a9 
+6(1 - 40fl'a"3 + a=^a(3)' 
(1-60' 



Pik,) 



-16a%" + 16oo'3o(=^) + (29 - 1080ao'2o"2 - Sa''a'''a^^^ 
2o3o"o(^) + 18(4^ - l)o2o'o"a(3) + 2o3o'o(^)] , 
96o'^o" - 96oo'=^o(3) - (209 - 5400ao''o"' + 48o'o"o(^) 



768A:27r2a9 

+ (34 - 120e)o2o"3 + 3(24^ - 7)0^0(3)' + 4(18^ - 7)o3o"o(^) + 2(89 - 2520a2o'o"o(3' 
-14oVo(^)+2o^o(6) 



(18) 



^ In [11], the stress-energy tensor for the massive field has been calculated in a covariant way from the 
quantum effective action obtained by point splitting, which must be equivalent to (17). 
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We check that both (17) and (18) obey the conservation equation p' + 3/i(p + P) = 0, as 
they should, since this is guaranteed by the adiabatic regularization. For each term in these 
expressions the total number of time derivatives acting on the scale factors always equals 
six and thus the magnitudes of p and P are fixed by jvn^ and H^a^ /k1, for massive and 
massless fields respectively, which is expected by dimensional analysis. Note that p{k^) and 
P{ktf) vanish identically for the conformally coupled scalar with ^ = 1/6. 



III. POWER LAW EXPANSION 

The final expressions (17) and (18) are not very illuminating. Therefore, in this section 
we focus on power law expansion in conformal time and set 

In terms of the proper time, which is defined as dt — adrj, the scale factor becomes 



a — 



-Y, a^^. (20) 



Since n = a/{a — 1), the metric linearly expanding in proper time with a = 1 must be 
analyzed separately, which we study at the end of this section. 



A. Massive case 

In the background (19), the stress-energy tensor of the massive field (17) becomes (recall 
that H — a' /a^ is the Hubble parameter in proper time) 

77—2 

p{m) = C„ P{m) = ^ p(m), (21) 



m?' ■ n 



where the constant C„ is given by 
1 



C — - 
" 201607r2n4 



3060 - 4742n - 402971^ + 471671^ + 1365n^ 
126(240 - 438n - 389n^ + 5n^(92 + 37n))^ + 3780(n + 1)(20 - 74n + 25n^ + 35n^)^^ 



-22680n(n + lf{lln - 10)C^J . (22) 

Prom the equation of state (21), we see that for < n < 2 the sign of the adiabatic vacuum 
pressure is the opposite of the energy density. This corresponds to the range a > 2 or a < 




FIG. 1. The graphs of n^Cn for ^ = (left) and = 1/6 (right). C„ is multiphed by since H 
contains a factor of n coming from a' . Although it cannot be seen in the graph for ^ = 0, the curve 
actually oscillates in the (—1,1) interval crossing the n-axis three times. 

in (20). Note that while for a > 2 the background is accelerating with decreasing Hubble 
parameter, for a < the Hubble parameter increases and there is a big-crunch singularity 
at finite proper time. 

It is interesting to compare the vacuum energy density with the background energy density 
driving the metric (19). From the Friedmann equation, we know that ~ pb/M^, where 
is the background energy density and Mp is the Planck mass. Thus one has 



Since adiabaticity requires m ^ H, unless H ^ Mp the vacuum energy density p(m) is 
much smaller than the background energy density ps showing that the backreaction effects 
can safely be ignored in this setup. 

For < 72 < 1, which corresponds to a < 0, the Hubble parameter H increases and 
thus adiabaticity will eventually be lost since H/m grows to become order one at some 
time before reaching the big-crunch singularity. Although the vacuum energy density p(m) 
also increases in time, when adiabaticity is lost, i.e. when H/m = 0{1), its ratio to the 
background energy density is given by H'^/Mp and thus p(m) is still negligible unless H is 
of the order of Planck scale. In any case, it would be interesting to study if the adiabatic 
vacuum can help to avoid the big-crunch singularity, along the lines of [12]. 

For minimally and conformally coupled scalars, i.e. for ,^ = and ^ = 1/6, the propor- 




(23) 



pB M'^m? ' 
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tionality constant C„ simplifies a lot: 

I 2oi607r2n4 ^060 - 4742n - 4029n^ + 4716n=^ + 1365n^] ^ = 0. 
In figure 1, we give plots of C„ for ^ = and ^ = 1/6. It can be seen from the graphs that 
there are intervals of the power n for which C„ is positive or negative. There also exists 
special values for which C„ vanishes. Therefore, depending on the expansion power, the 
(leading order contribution to) the vacuum energy density can be positive, negative or even 
zero. 

Let us finally focus on important special cases. In de Sitter space, which corresponds to 
n = 1, the equation of state becomes P{m) = —p{m) and thus adiabatic vacuum energy 
density is equivalent to a cosmological constant. In this case C„ is given as 

_ 270 - 7308^ + 45360f - 90720^^ 

- 20160;^^ ■ ^^^^ 

For ^ = 0, Ci > and for ^ = 1/6, Ci < 0, thus while the vacuum energy density is positive 
for the minimally coupled scalar, it is negative for the conformally coupled case, modifying 
the background cosmological constant correspondingly. It is interesting to compare these 
findings with the massive scalar field in Bunch-Davics vacuum studied in [14]. For m ^ H, 
the vacuum energy density in Bunch-Davies vacuum is given by pbd — —rn^ (see [13] or eq. 
(3.15) of [14]). Therefore, the vacuum energy densities corresponding to Bunch-Davies and 
adiabatic vacua differ both in magnitude and in sign. 

In the radiation dominated background with n = —1, the vacuum pressure satisfies 
P{m) = 3p(m), which corresponds to a very stiff matter. In that case the constant Cn 
becomes 

_ 211 - 882e 

- 100 807r2 ■ ^^^^ 
The vacuum energy density in the radiation dominated universe turns out to be positive 

both for the minimally and the conformally coupled scalars. 

Finally, in the matter dominated background one has n = —2 and P{m) = 2p(m), again 

equivalent to a very stiff matter. The constant C„ becomes 

36(3e-l)(29 + 960e-139 

- looso;^^ ■ ^^^^ 

One can see that both for ^ = and for ^ = 1/6, C_2 < and the vacuum energy density 
becomes negative. 
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B. Massless case 



In the background (19), the stress-energy tensor (18) corresponding to the UV modes of 
the massless field with k > becomes 

H^a^ n — 6 

p{h) = Dn P{h) = p{k.), (28) 

where the constant D„ is given by 

1 



5(l-60V(n+l)(2-n)(2 + n(n+l)(6^-l)) . (29) 



In that case the adiabatic vacuum pressure has the opposite sign compared to the vacuum 
energy density for < n < 6, which is equivalent to a < or a > 6/5. On the other hand, 
for < n < 3/2, which corresponds to ck > 3 or a < 0, the vacuum energy density grows 
and adiabaticity will be lost in time since Ha/k^ grows. Interestingly, in de Sitter space 
with n = 1, while the stress-energy tensor vanishes identically for the conformally and the 
minimally coupled scalars, the energy density grows for a generic ^, i.e. it is not in the form 
of a cosmological constant. In that case however, the vacuum energy density still stays much 
smaller than the background energy density as long as the adiabatic approximation holds . 

It is important to emphasize that for the high energy modes with k ^ aif, the adiabatic 
vacuum is equivalent to Bunch-Davies vacuum, which has the following mode functions 

P^k^'-V^hM, (30) 
where hu denotes spherical Hankel function of first kind and 



u^-l + ^l + n(n + l)(l-60. (31) 

Using the asymptotic form of the spherical Hankel functions, one can find as A; — > oo that 

„BD.t^,-^kr,(. + ^) uju + 1) [{1 - GQnjn + 1) - 2] \ 



which exactly coincides up to a phase with the adiabatic expansion (this can easily be checked 
for the the two terms written above). Note that for some special values, hke n — —1, the 
series terminate and one ends up with a finite series instead of getting infinitely many terms. 
On the other hand, for infrared (IR) modes adiabatic expansion fails but (30) can still be 

11 




FIG. 2. The graph of n^Dn for ^ = 0. D„ is multipUed by since H contains a factor of n coming 
from a'. 



used for Bunch-Davies vacuum (however the vacuum expectation values are now plagued by 
IR divergences, see [7]). 

For the minimally coupled scalar, the Dn coefficient becomes 

5(n^ -Sn^ + 4) 



Dr. 



1287r2n4 



(33) 



Curiously, the stress-energy tensor vanishes in the de Sitter space (n = 1) and in the radiation 
dominated (n = —1) and the matter dominated {n = —2) universes (it also vanishes for n = 2 
corresponding to a{t) = in proper time). As before, depending on n the vacuum energy 
density can be positive or negative, which can be seen from the plot of Dn given in figure 2. 



C. The linear expansion a = 1 



In the linearly expanding universe where a{t) = t/to, the scale factor in the conformal 
time becomes 



a = exp 



^0 



(34) 



Using (17), one can see that the stress-energy tensor of the massive scalar obeys P{m) = 
p{m) where 



On the other hand, the stress-energy tensor for the massless UV modes satisfies P{k^: 
p(A;*)/3, where 



5(1-60^ 
128A;27r2 



(36) 
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In both cases, the stress-energy tensor vanishes for the conformally coupled scalar and it 
is positive for the minimally coupled one. Note that the equation of state can be obtained 
from (21) and (28) by taking n — >■ oo limit, which corresponds to a = 1. 

IV. CONCLUSIONS 

In this paper, we determine the leading order contribution of the adiabatic mode func- 
tions to the vacuum expectation value of the stress-energy tensor by explicitly calculating 
adiabatic order six terms. Due to the assumption of adiabaticity, the particle creation ef- 
fects are absent in this setup and thus the resulting expressions can be thought to be related 
to vacuum polarization effects. When the scale factor of the universe is given by a simple 
power, the stress-energy tensor simplifies a lot and it obeys a simple equation of state. In 
such a background, the vacuum energy density is proportional to a single term, which can 
actually be fixed by dimensional analysis. We determine when it is possible to use adiabatic 
approximation and when it fails in time. Depending on the power of the scale factor, (the 
leading order contribution to) the vacuum energy density can be positive, negative or zero. 
Moreover, it can be observed to increase, decrease or stay constant for different values of 
the power. 

For a massive field it is natural to assume adiabatic vacuum when H. For example 
at early times, the stabilized moduli, which are known to exist in string theory motivated 
models, must satisfy this condition (this is usually required for not to change the standard 
cosmological picture, which is known to be in good agreement with observations). At late 
times, nearly all known massive fields (neutrinos can be an exception) obey this condition and 
thus assuming adiabatic vacuum for them is reasonable. For massless fields, only the high 
energy modes can be thought to evolve adiabatically, where the critical scale is determined 
by the Hubble parameter. Moreover, the adiabatic vacuum coincides with the Bunch-Davies 
vacuum for UV modes, at least when the scale factor is a simple power (although we are 
not aware of any proof, this is plausibly true for any given scale factor). 

We see that in the adiabatic vacuum, the energy densities corresponding to a very massive 
scalar field and the UV modes of a massless scalar field turn out to be very small compared 
to the background energy density driving the expansion (As noted earlier that this is true 
even when the leading order contribution to the vacuum energy density grows in time and 
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adiabaticity is eventually lost. However, it would be interesting to study the evolution of the 
vacuum energy density in such a case.) This supports the idea that the cosmological constant 
problem is an IR issue rather than being a UV problem, which is contrary to the common 
thought emphasizing the huge order of magnitude discrepancy between the expected and 
observed values of the cosmological constant. Namely, with a proper regularization the 
contributions of the high energy modes to the vacuum energy density become negligible. 
Physically the situation is very similar to the Casimir energy associated with parallel plates 
where the vacuum energy density turns out to be fixed by the distance between plates, i.e. 
by the IR scale rather than by the UV scale (moreover it turns out for Casimir energy that 
the massive fields play no role and can safely be ignored). Since experimentally measured 
Casimir energy is in good agreement with theoretical calculations, the previous statement 
can be claimed to have a sound basis. Of course, the IR problem can be seen to be more 
challenging, but in any case it is important to pin down the real issue and our findings 
support the idea that IR physics plays the key role in the cosmological constant problem. 
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